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Abstract 

The goal of this paper to determine the optimal horoball packing ar- 
rangements and their densities for all four fully asymptotic Coxeter tilings 
(Coxeter honeycombs) in hyperbolic 3-space H 3 . Centers of horoballs are 
required to lie at vertices of the regular polyhedral cells constituting the 
tiling. We allow horoballs of different types at the various vertices. Our re- 
sults are derived through a generalization of the projective methodology for 
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hyperbolic spaces. The main result states that the known Boroczky-Florian 
density upper bound for "congruent horoball" packings of H 3 remains valid 
for the class of fully asymptotic Coxeter tilings, even if packing conditions 
are relaxed by allowing for horoballs of different types under prescribed 
symmetry groups. The consequences of this remarkable result are discussed 
for various Coxeter tilings. 



1 Introduction 

Local optimal ball packings for regular tilings have been studied extensively in 
the literature. Of special interest are tilings in hyperbolic n-space; Boroczky and 
Florian HB-F64M gave the universal density upper bound for all congruent ball 
packings of the 3-dimensional hyperbolic space H 3 without any symmetry as- 
sumptions. This classical result provides the density upper bound realized by a 
regular horoball packing of (3, 3, 6) in H 3 , as shown in Section 2. The optimal 
density is related to the Dirichlet-Voronoi cell of every ball, as follows: 

«o = (l + i-^-p + i + p- - + +...)" 1 * 0.85327609. 

This limit is achieved by 4 horoballs centered at the vertices of a regular ideal sim- 
plex, tangent to each other at the "midpoints" of the edges, i.e. as the projection 
of the simplex center into any edge of the simplex. 

Papers del, HBo-RL IKj-K-KI . flEED, flRMTl introduce recent novel develop- 
ments in the classical topic of the ball (or sphere) packings of H 3 . Locally dense 
(optimal) ball, horoball packings in H 3 are of great significance, as important in- 
formation regarding crystal structures can be obtained using locally optimal ball 
and horoball arrangements. 

The present work is based on the projective interpretation of the hyperbolic ge- 
ometry, proposed in [M93 1, [M97 1. In subsequent works the second author studied 
a class of face transitive tilings HD-H-ML the so-called generalized Lambert-cube 
tilings in HSz03-l|| and HSz05-li where an algorithmic approach for determining 
volumes of hyperbolic polyhedra was developed and implemented. Using this 
novel approach, the locally optimal ball packings were found for the configura- 
tions in which the ball and horoball centers lie either within the Lambert-cubes 
or at the vertices of the cubes, respectively, and the optimal packing densities 
of the corresponding tilings were computed [Sz 05-li Optimal ball and horoball 
packings of the regular Coxeter honeycombs in M d , (d > 3) with one horoball 
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type were found in HSz05-2ll . HSz07-lH . and the optimal ball and horoball packing 
densities computed. 

In related work, t/-dimensional (d > 3) hyperbolic prism honeycombs gen- 
erated by "inscribed hyperspheres" were investigated in HSz06-ll and [Sz06-2J. 
The optimal hyperball packings of infinitely many 3-dimensional prism tilings 
(mosaics) together with their metric data were determined in HSz06-21l . In hyper- 
bolic 4-space H 4 there are only 2 honeycombs with metric data corresponding to 
their 3-dimensional counterparts. The densities of the optimal hyperball packings 
in 4-space are determined in HSz06-2H . In H 5 there are 3 types of such mosaics, 
and the corresponding problems are extensively studied. In M d (d > 5) there are 
no longer any regular prism tilings. 

In this paper, we study locally optimal ball and horoball packings in the four 
fully asymptotic hyperbolic tilings of H 3 , while allowing different types of horoballs 
to be centered at the vertices of the honeycombs. In Section 2, we provide pre- 
liminaries on the rf-dimensional honeycombs. In Section 3, we introduce the pro- 
jective model [Sz05-1 1 to determine the densities of the optimally dense horoball 
packings in hyperbolic space M. d . In Section 4 we determine the optimal packing 
densities in H 3 for various honeycombs, when horoballs of various types are al- 
lowed. We find that the densest possible packings yield density values identical 
to that of the Boroczky-Florian bound HB-F641 . In all studied configurations, the 
optimal densities never surpass the Boroczky-Florian upper bound, even when 
replacing the "congruency" constrains with regularity constraints. We finish the 
paper with conclusions and directions for future research. 



2 Overview on d-dimensional hyperbolic 
honeycombs 

Hyperbolic geometry is based on the principles of Bolyai-Lobachevsky geometry 
[P06J. A c/-dimensional honeycomb V, also referred to as a solid tessellation or 
tiling, is an infinite collection of congruent polyhedra (polytopes) that fit together 
face-to-face to fill the entire geometric space (H d (d ^ 2)) exactly once. We take 
the cells to be congruent regular polyhedra. A honeycomb with cells congruent 
to a given regular polyhedron P exists if and only if the dihedral angle of P is 
a submultiple of 2n (in the hyperbolic plane zero angles are also permissible). 
A complete classification of honeycombs with bounded cells was first given by 
SCHLEGEL in 1883. The classification was completed by including the polyhedra 
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with unbounded cells, namely the fully asymptotic ones by Coxeter in 1954 
HC56H . Such honeycombs exist only for d < 5 in hyperbolic <i-space M d . In this 
paper Coxeter honeycombs or Coxeter tilings refer to tilings described in Table 1. 

An alternative approach to describing honeycombs involves analysis of their 
symmetry groups. If V is a Coxeter honeycomb, then any rigid motion moving 
one cell into another maps the entire honeycomb onto itself. The symmetry group 
of a honeycomb is denoted by SymV. The characteristic simplex T of any cell 
P E V is a fundamental domain of the symmetry group SymV generated by 
reflections in its facets which are (d — 1) -dimensional hyperfaces. 

The scheme of a regular polytope Pisa weighted graph (diagram) characteriz- 
ing P C tf up to congruence. The nodes of the scheme, numbered by 0,1, ... ,d, 
correspond to the bounding hyperplanes of J-'. Two nodes are joined by an edge if 
the corresponding hyperplanes are non-orthogonal. Let the set of weights (m, n 2 , 
ris, . . . , rid-i) be the Schlafli symbol of P, and rid be the weight describing the 
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dihedral angle of P, such that the dihedral angle is equal to — . In this case T is 
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Figure 1: Coxeter-Schlafli simplex scheme 

The Schlafli symbol of the honeycomb V is the ordered set (ni, n 2 , n 3 , . . . , 
n d-i, n d) above. A (d + 1) x (d+ 1) symmetric matrix (6 ,J ) is constructed for 
each scheme in the following manner: b u = 1 and if i ^ j E {0,1,2, ... ,d} then 
W = — cos — . For all angles between the facets i,j of J 7 holds then n^ = n^-i k- 

riij 

Reversing the numbering of the nodes of scheme V while keeping the weights, 
leads to the scheme of the dual honeycomb V* whose symmetry group coincides 
with SymV. 

In this paper we investigate regular Coxeter honeycombs and their optimal 
horoball packings in the hyperbolic space H 3 , where the horoballs are allowed 
to be of different types. SymV denotes the symmetry group of the honeycomb 
V thus 

J^n 1 n2...n d — \ {_) 7Vr n\n,2—n,d) }■ 

7 e SymV ni „ 2 ...n d _ 1 

In order to calculate the packing density, we relate each ball or horoball, respec- 
tively, to its regular polytope P nin2 ... nd in which it is contained. These polytopes 
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are not necessarily assumed to be Dirichlet-Voronoi cells. 

Table 1: Classification of 3 -dimensional Coxeter tilings 



No. 


Description 


Schlafli symbol (p, q, r) 


1. 


Cells having proper centers and vertices 


(3,5,3), (4,3,5), (5,3,4), (5,3,5) 


2. 


Fully asymptotic cells 


(3,3,6), (3,4,4), (4,3,6), (5,3,6) 


3. 


Infinite centers and proper 
or ideal vertices 


(3,6,3), (4,4,4), (6,3,6) 
(4,4,3), (6,3,3), (6,3,4), (6,3,5) 



As listed in Table 1, Coxeter tilings with parameters in row 1 include cells 
having proper centers and vertices. The polyhedra of honeycombs with Schlafli 
symbols in row 3 of Table 1 , have infinite centers and proper or ideal vertices. The 
polyhedra of tilings in row 2 of Table 1 are called fully asymptotic; moreover their 
centers are proper and their vertices lie on the absolute of the hyperbolic space i.e. 
they are ideal vertices. 



3 The projective model 

3.1 Basic Notions 

Let X denote one of either the d-dimensional sphere E> d , the d dimensional Eu- 
clidean space E d , or the hyperbolic space H d , d > 2. For M d we use the projective 
model in Lorentz space E 1 ^ of signature (1, d), i.e. E 1,(i is the real vector space 
yrf+i e q U ipp e d vvith the bilinear form of signature (1, d) 

(x, y) = -x°y° + x^ 1 + ■ ■ ■ + x d y d (3.1) 

where the non-zero vectors 

x = (x°, x\...,x d )e V d+1 and y = (y°, y 1 ,..., y d ) e V d+1 , 

are determined up to real factors and they represent points in V d (K). M d is repre- 
sented as the interior of the absolute quadratic form 



Q = {[x] GP d |(x, x) = 0} = <9H fl 



(3.2) 



in real projective space V d (Y d+l , V^+i). All proper interior point x G M d are 
characterized by ( x, x) < 0. 
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The points on the boundary dM. d in V d represent the absolute points at infinity 
of M d . Points y with ( y, y) > lie outside dM d and are called outer points of 
M d . Let -P([x]) G V d ; a point [y] G V d is said to be conjugate to [x] relative to 
Q when ( x, y) = 0. The set of all points conjugate to P([x]) form a projective 
(polar) hyperplane 

pol(P):={[y]eV d \(x, y) = 0}. (3.3) 

Hence the bilinear form Q by (3.1) induces a bijection (linear polarity V d+1 — > 
Vd+i)) from the points of V d onto its hyperplanes. 

Point X[x] and the hyperplane a [a] are called incident if the value of the 
linear form a on the vector x is equal to zero; i.e., xa = (x G V d+1 \ 
{0}, a G V d+ i \ {0}). Straight lines in V d are characterized by the 2-subspaces 
of V d+1 or (d - l)-spaces of V d+1 HM97H . 

Let P C M d denote a polyhedron bounded by hyperplanes H\ which are 
characterized by unit normal vectors b l G Vd+i directed inwards with respect to 
P: 

W := {x G M d \( x, V) = 0} with (b\ V) = 1. (3.4) 

We always assume P to be an acute-angled polyhedron and the vertices to be 
proper points or to lie at infinity. 

The Gram matrix G{P) := {{b\ b j )) i,j G {0, 1, 2 . . .d} of normal vectors V 
associated with P is an indecomposable symmetric matrix of signature (1, d) with 
entries (b\b l ) = 1 and (b l ,b>) < for i ^ j, having the following geometrical 
meaning 



(b\V) 



(0 if H' l ±H j , 

— cos q' j if H l , H 3 intersect on P at angle a %3 , 

— 1 if H l , H 3 are parallel in hyperbolic sense, 

— cosh f J if H l , H 3 admit a common perpendicular of length l %3 . 



Definition 3.1 [K911 IB-Hl An orthoscheme O in X is a simplex bounded by d+ 1 
hyperplanes H°, . . . , H d such that 

W-LH 3 , forj^i-l,i,i + l. 

A plane orthoscheme is a right-angled triangle, the area of which can be ex- 
pressed by the defect formula. For orthoschemes we denote the (d — l)-hyperface 
opposite to the vertex Ai by H % (0 < i < d). An orthoscheme O has d dihedral 



Packings for Fully Asymptotic Coxeter Tilings. . . 7 

angles different from right angles. Let a u denote the dihedral angle of O between 
the faces H % and ffl. Then we have 



Ti 



a' 



-, if 0<i< j-1 <d. 



The remaining d dihedral angles cn v+1 , (0 < i < d—1) are called the essential 
angles of O. The initial vertex A and final vertex A d of the orthogonal edge-path 



d-l 



U AiAi+t 



i=0 



are called principal vertices of the orthoscheme. 

In this work, the characteristic simplex T of any honeycomb V with Schlafli 
symbol (ni,n 2 , n 3 , . . . , n d ) is an orthoscheme. 

The matrix (6* J ) = G(P) is the so called Coxeter-Schlafli matrix of the or- 
thoscheme T with parameters n 1 , n 2 , ra 3 , . . . , n d : 



(b ij ) :-■ 



/ 1 

— cos 





»i 



— cos 

1 

— cos 





"1 



ll 2 





— cos 

1 

— cos 



"2 



»:i 



o\ 







V 



o 



o 



— cos 



(3.5) 



rid / 



Inverting the Coxeter-Schlafli matrix (W) (3.5) of an orthoscheme gives the 
matrix (oy), which can be used to express distances between two vertices through 
the formula [Sz06-2]: 



cosh — — 
k 



v 7 



"a jj 



In this paper we set the sectional curvature of M d , K = —k 2 , to be k = 
distance s of two proper points (x) and (y) is calculated by the formula: 



cosh 



-<x, y) 



k v/(x, x)(y, y) 



(3.6) 
1. The 

(3.7) 



8 R. T. Kozma and J. Szirmai 

3.2 Characterization of horoballs in Hyperbolic Space H 3 

Definition 3.2 A horoshere in the hyperbolic geometry is the surface orthogonal 
to the set of parallel lines, passing through the same point on the absolute quadratic 
surface (simply absolute) of the hyperbolic space. 

We represent hyperbolic space H 3 in the Cayley-Klein ball model. We intro- 
duce a projective coordinate system using vector basis bj (i = 0, 1, 2, 3) for V 3 
where the coordinates of center of the model is A 2 (l, 0, 0, 0). We pick an arbitrary 
point at infinity to be A 3 (l, 0, 0, 1). 

As it is known, the equation of a horosphere with center A 3 (l, 0, 0, 1) through 
point S(l, 0, 0, s) is derived using the surface pencil of the absolute sphere and a 
plane tangent to the sphere at point A 3 (l, 0, 0, 1). The equation of the absolute 
sphere is — x°x° + x 1 x 1 + x 2 x 2 + x 3 x 3 = 0. The equation of a plane tangent to 
the absolute of our model at point A 3 (l, 0, 0, 1) is x° — x 3 = 0. 

The general equation of the horosphere is 

= A(-:r :r + x 1 x 1 + x 2 x 2 + x 3 x 3 ) + /2(x° - x 3 f . 
This passes through point S(l, 0, 0, s) so we may write 

A(-i + s 2 )+M-i + *) 2 = o^-- (s lr 



jl 1 — s 2 

If s ^ 1, then {S ~ { (-x°x° + x l x l + x 2 x 2 + x 3 x 3 ) + (x° - x 3 f = <£> 

1 — s z 

&(s- l)(-x°x° + x l x l + x 2 x 2 + x 3 x 3 ) - (1 + s)(x° - x 3 ) 2 = 

This way we obtain the following equation for the horosphere in our Cayley- 
Klein model of H 3 : 

- 2sx°x° - 2x 3 x 3 + 2(s + l)(x°x 3 ) + (s - l)(:rV + x 2 x 2 ) = (3.8) 



Remark 3.3 We have obtained the equation of the horosphere in the Cartesian 
coordinate system: (x := 4, y := 4, z := 4) 



2(x 2 + y 2 ) 4(s-^) 
l-s + (l-sf 



(3.9) 
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Remark 3.4 It is useful for visualization purposes to convert the horosphere equa- 
tion into polar coordinates. By multiplying the polar coordinate form by rotation 
matrices we can easily obtain horospheres around arbitrary points at infinity in the 
model. The polar form is by parameters <\> E [0, 2ir), 9 E [0, n] 



x — \l sin 9 cos y = \ sin 9 sin < 



2 



1+s 1-s 
z = 1 cos 9. 



(3.10) 



3.3 Volumes of horoball sectors 

The length l(x) of a horospheric arc of a chord segment x is determined by the 
classical formula due to J. BOLYAl: 

x 
l(x) = /csinh— (at present A; = 1). (3.11) 

k 

The intrinsic geometry of the horosphere is Euclidean, therefore, the area A of 
a horospherical triangle is computed by the formula of Heron. The volume of the 
horoball pieces can be calculated using another formula by J. Bolyai. If the area 
of a domain on the horoshere is A, the volume determined by A and the aggregate 
of axes drawn from A is equal to 

V = -kA (we assume that k — 1 in this paper). (3.12) 

4 Horoball packings for totally asymptotic Coxeter 
honeycombs 

4.1 Basic results on packing in H 3 

In this section we determine optimal horoball packings for the four totally asymp- 
totic Coxeter tilings V pqr with Schalfi symbols (p, q, r) = (3, 3, 6), (3, 4, 4), (4, 3, 
6), (5, 3, 6). Vertices of a regular cell P pqr are denoted by Et, (i — 0, 1, 2, 3, 4 . . . ), 
and lie on the absolute of H 3 , hence these vertices are centers of horoballs. The 
number of the vertices of P pqr is denoted by N pqr and we write Bi for the horoball 
centered at E^ . We require a horoball to lie at every ideal vertex of the honey- 
comb, and we vary the touching types of the horoballs. 
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Remark 4.1 For example, if (j>, q, r) = (3, 3, 6) then we obtain the "tetrahedral 
case" where N 33e = 4. For parameter (p, q, r) = (3, 4, 4) we get the "octahedral 
case" and N a u = 6. 

The type of a horoball is allowed to expand until either the horoball comes 
into contact with other horoballs or a non-adjacent side of the honeycomb. These 
conditions are satisfactory to ensure that the balls form a non-overlapping horoball 
arrangement, as such the collection of all horoballs is a well defined packing in 
H 3 , denoted by B pqr . 

Definition 4.2 The density of a horoball packing in Coxeter honeycomb V pqr is 
defined as 

6{Bmr) = SST^niW) 



Vol(P pqr ) 

The aim of this section is to determine the optimal packing densities for the 
four totally asymptotic tilings (see Table 1) in 3-dimensional hyperbolic space H 3 . 
We will make heavy use of the following Lemma [Sz 05-lH : 

Lemma 4.3 Let there B\(x) and B 2 (x) denote two horoballs with centers C\ and 
C 2 respectively. Take T\ and r 2 to be two congruent trihedra, with vertices C\ 
and C 2 . Assume these horoballs B\ and B 2 are tangent at point I(x) G C1C2 
and C1C2 to be a common edge of the two trihedra T\ and r 2 . We define the point 
of contact 1(0) such that the following equality holds for volumes of horoball 
sectors: 

v(o) := 2Voi(B l (o) n n) = 2^0/(52(0) n r 2 ). 

Ifx denotes the hyperbolic distance between 1(0) and I(x), then the function 

V(x):= Vol(B x (x) n n) + Vol(B 2 (x) n r 2 ) 

strictly increases as x — > ±00. 

Proof: Let £ and C be parallel horocycles with centre C and let A and B be 
two points on the curve C and A' := CA fl £', B' := CB fl C. By the classical 
formula of J. BOLYAI 

n(A'B') 

H(AB) 
where the horocyclic distance between A and B is denoted by %(A, B). 
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Figure 2: C\ and C 2 are centers of the horoballs B\ and B 2 , and I(x) is the point 
of contact. 1(0) represents the point at which Vol(Bi) = Vol(B 2 ). 



Then by the above formulas we obtain the following volume function: 
V(x) = Vol(B x {x) n n) + Vol(B 2 (x) n r 2 ) = 
= iv(0)(e* + -i)=y(0)ooah(^). 

It is well known that this function strictly increases in the interval (0, oo). This 
Lemma is illustrated in Fig. 2. □ 

Corollary 4.4 The statement of the Lemma holds for horoballs intersecting the 
arbitrary two congruent frames of rays joining C\ and C 2 above, respectively. 

Proof: Follows from the Lemma 4.3, by dividing the frames of rays into con- 
gruent trihedra. □ 



4.2 The (3, 3, 6) tetrahedral tiling 

The (3, 3, 6) Coxeter tiling is a three dimensional honeycomb with cells comprised 
of fully asymptotic regular tetrahedra. We arbitrarily select one such tetrahedron 
E EiE 2 E 3 (see Fig. 3), and place the horoball centers at vertices E , . . . , E 3 . We 
vary the types of the horoballs so that they satisfy our constraints of non-overlap. 
The packing density is obtained by Definition 4.2. 

Define the orthoscheme A AiA 2 A 3 as follows: A = E and A 3 = E 3 are 
two vertices of the tetrahedra (see Fig. 3); A 2 is the center of the triangle E E 1 E 2 
opposite the vertex A 3 , and A\ is the footpoint of A 3 on the edge EqE\. One 
tetrahedral cell is decomposable into 6 such congruent orthoschemes. The Schlafli 
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symbol of ortho scheme A AiA 2 A 3 is (3, 6, 3), and the orthoscheme is labeled by 

O (3,6,3)- 

B and B 3 are two horoballs centered at E and E 3 , i.e., the two vertices of 
the tetrahedra common with the orthoscheme. The density of the (3, 3, 6) Coxeter 
tiling is obtained using Definition 4.2: 

Xftt ^ _ V0l{B Q n (9(3,6,3)) + V0l(B 3 H (9(3,6,3)) 
^1^336 J — 



^(0(3,6,3): 



Proposition 4.5 The packing density obtained in 0(3,6,3) can be extended to tetra- 
hedron P336 and therefore to the entire H 3 . 

Proof: We consider the following steps: 

1 . In an optimally dense packing, at least two horoballs must touch each other 
in the tetrahedron, otherwise the density could be improved by blowing up 
any one horoball until it touches a neighboring horoball. 

2. If two horoballs touch at the "midpoint" of edge A A 3 as projection of the 
simplex centre on it, then by blowing up the remaining two horoballs, they 
will also touch at "midpoints", due to symmetry considerations. Note, that 
this case is the arrangement in the Boroczky-Florian density upper bound 
using the same horoballs. 

3 . Given two horoballs tangent at a non-midpoint of an edge, then the horoball, 
having the midpoint in its interior, will contain the "midpoint" of all 3 edges 
extending from its center. As a result the remaining 3 horoballs should be 
of the same touching type. The "small horoballs" each is tangent only to 
the large horoball. 

We just showed that the "largest horoball" determines the configuration of all 
other horoballs, and as a consequence the packing density. One parameter corre- 
sponding to the "largest horoball" suffices to determine the packing density for all 
candidates of optimal density. 

Due to symmetry, it is enough to consider cases within orthoscheme (9( 3i6)3 ), 
where the horoball at E expands from the midpoint until it becomes tangent to 
the side of the cell opposite to it. Assume no balls cover the midpoint along 
E E 3 . Then the horoball B can be expanded until the midpoint. In this case, 
the horoball B 3 is contained within 0(3 j6;3 ). Finally, the packing density can be 
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varied by expanding the horoball B while keeping B 3 tangent to it. If we expand 
Bq until it touches A 2 , all candidates for optimal packings within the orthoscheme 
are considered. The densities obtained from the orthoscheme 0(3,6,3) will cover 
all candidates for optimally dense packings of honeycomb (3,3,6). Densities 
determined within the orthoschemes can be generalized to the entire H 3 by the 
symmetries of P336. □ 

In the rest of this section, we prove the basic theorem on the optimal packing 
density in (3, 3, 6). First, we define the tangent point 1(0) E A A 3 of horoballs 
B and B 3 so that the following equality holds for the volumes of the horoball 
sectors 

V(0) := 2Vol(B (0) n (3 , 6>3 )) = 2Vol(B 3 (0) n 0(3,6,3))- 

Consider point I(x) on the edge A A 3 where the horoballs B^x), (i = 0, 3) 
are tangent at point I(x) E A A 3 . as the hyperbolic distance between 1(0) and 
I(x). Function V(x) is defined as follows: 

V(x) := Vol(B (x) n (9(3,6,3)) + Vol(B 3 (x) H 0(3,6,3))- 

By Lemma 4.3 it follows that function V(x) strictly increases as I(x) (x E 
[— arctanh(l/2), arctanh(l/2)]) moves away from 1(0) along A A 3 . That implies 
that the density function 



xf „ w, V(x) Vol(B (0)nO {3A3) )e 2x + Vol(B 3 (0) nO (3 , 6 , 3 ))e 



-2x 



0(3,6,3) V"o/(0(3, 6 , 3 )) 

attains its maximum at the two endpoints of the interval [— arctanh(l/2), arctanh(l/2)]. 

Definition 4.6 A coordinate system is assigned to the orthoscheme A AiA 2 A 3 ; 
let A 2 := (1, 0, 0, 0) be the origin, A := (1, 0, 1, 0), A 3 := (1, 0, 0, 1) and A x = 
(1, ^, i, 0) where A x is the "midpoint" of the edge E E 1 (see Fig. 3). 

In the orthoscheme AqAiA 2 A 3 , the horoball centers are located at points A 
and A 3 , and the horoballs meet along side A A 3 (See Fig. 3 (i)). In order to 
determine the properties of the packing, we calculate the five points of intersection 
of the horoballs with the edges of the orthoscheme. The distances of these five 
points of intersection determine the area of the horospheric triangles, hence the 
volume of the horoball sectors through Bolyai's formulas (3.11-12). 



14 



R. T. Kozma and J. Szirmai 




Figure 3: a. Orthoscheme 0(3,6,3); b. the totally asymptotic regular tetrahedron 
E E 1 E 2 E 3 for(3,3,6). 



In order to simplify the computations and visualizations of horospheres we 
introduce a new parameter s. We define s to be the Euclidean parameter of the 
origin A 2 and point 5(1,0,0, s) along segment A 2 A 3 . Based on the proof of 
Proposition 4.4, s is the only parameter that determines the density of the packing. 

All horoball configurations that yield valid candidates for the optimal packing 
density occur while we continuously vary the horoball parameter s G [0, 1/2]. In 
order to determine the optimal packing, we proceed as follows. We express the 
density 5(B3sq) as the function of s, study its behavior and determine its extremal 
points. If s = 1/2, all horoballs are in same type, thus we achieve the packing 
arrangement B\ m , which is the Boroczky-Florian case with the known packing 
density. On the other hand, for s = 0, we find a different horoball arrangement 
i3| 36 (see Fig. 3) with the same density as for s = 1/2. Finally, from Lemma 4.3 
it follows that all other densities with horoball parameter s G (0, 1/2) are smaller. 

First we calculate the five intersections of the edges of 0(3,6,3) and the two 
horoballs B and B 3 for s = 0. Recall that these are all a function of the "type" of 
the large horoball, and depend on the parameter s. The length of the sides of the 
horospherical triangle, with vertices Xi, Yi, Z% £ Bi, i E {0, 3} are given by 



2sinh 



d(Xi, Yi t 



hi := 2sinh 



d{Yi,Zi) 



2sinh 



d{Zi,Xi) 



Based on the side lengths, Heron's formula can be used to obtain the area 
Ai(s) of horospherical triangles XjYjZj. Using Bolyai's volume formula (3.12) 
for horoball pieces, the Definition 4.2, and setting k — 1, the density of the pack- 
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0.0 0.1 0.2 0.3 0.4 0.5 



a. s E (0, 1/2) b. s = 1/2 c. s = 

Figure 4: Horoball packings of (3, 3, 6) tiling. 



ing can be expressed as a function of s as well 



5(B 336 )(s) 



l(A (s) + A 3 (s)) 
Vol{O m3) ) 



se [0,1/2]. 



From Lemma 4.3 it follows, that the optimal densities are realized on the end- 
points of the interval [0, 1/2]. In order to determine the highest packing density, 
we calculate the density for horoball pieces with s — 0. The volumes of or- 
thoscheme 0(3,6,3) in H 3 are calculated using Lobachevsky's volume formula (see 
HSzQ3-llO . All other densities with < s < 1/2 can be evaluated using the corre- 
sponding volume formula. The results are displayed in Fig. 4. a, while Fig. 4.b and 



Fig. 4.c illustrate the ball arrangements B 



\ m (i = 1, 2) for s = 1/2 and s = 0, 
respectively. Note that B\ m and B 336 are two different cases. In the first case the 
horoballs are of same type and they touch each other pairwise, while for s = we 
have balls with 2 distinct types, and the "small balls" do not meet. As a result, we 
have just proved: 

Theorem 4.7 There are two distinct optimally dense horoball arrangements B l 336 , 
(i — 1, 2) for the tetrahedral Coxeter tiling (3, 3, 6) with the same density: 5(B l 336 ) : 
0.85327609. 



4.3 The (3, 4, 4) octahedral tiling 

In this section we consider horoball packings with centers located at ideal vertices 
of the octahedral honeycomb (3,4,4). Our approach is as in the previous sec- 
tion for the (3, 3, 6) Coxeter tiling. We again allow the horoballs to be of different 
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types. The (3, 4, 4) Coxeter tiling decomposes the 3-dimensional hyperbolic space 
into congruent cells consisting of regular fully asymptotic octahedra. Four octahe- 
dra meet along each edge. As in the case of (3, 3, 6), we again choose one cell of 
the tiling in order to perform our density calculations. Again, we vary the types of 
the horoballs centered in the honeycomb vertices. The packing density obtained 
form this cell is again extended to H 3 by SymV 3 ^. 

We define the orthoscheme AqA\A 2 A 3 of an octahedron for the calculations. 
Let A and A 3 be two adjacent vertices of the octahedron E EiE 2 E 3 E4E 5 (see 
Fig. 5. a-b-c), A 2 be the center of the octahedra, and take A\ as the midpoint, in 
the Euclidean sense, of the edge extending from A sharing a common facet with 
A3. This orthoscheme is C444 and has Schlafii symbol (4, 4, 4). 

Definition 4.8 The following coordinate system is defined for orthoscheme 
A AxA 2 A 3 (seeFig. 6): 

A 2 = (1, 0, 0, 0), A = (1, 0, 1, 0), A 3 = (1, 0, 0, 1), Ax = (1, ^ ^ 0). 

Using Lobachevsky volume formula for orthoschemes, we obtain the volume 
of one octahedron HSz05-2ll : Vol(P U4 ) = 16 • Vol(0 {4AA) ) « 3.66384. Applying 
the definition of the packing density for the case of tiling (3, 4, 4), we obtain: 

HBui) = azfnM , ( 4.i) 

Vol(P 3 u) 

where B,i fi -P344 (i — 1, . . . , 6) denote the 6 horoball sectors, one in each vertex of 
the octahedron P344, and we assume that the horoballs Bi form a horoball packing 
in H 3 . 

We consider the following three basic horoball configurations £>g 44 , (i = 1,2,3): 

1. All 6 horoballs are of the same type and the adjacent horoballs touch each 
other at the "midpoints" of each edge. We define the point of tangency of 
two horoballs B and B 3 on side A A 3 to be 1(0) so that the following 
equality holds: 

7(0) := 6 ■ Vol(B (0) H P344)) = 6 • Vol(B 3 (0) H P 344 ) = 6 ■ V . 

In this case V := Vol(Bi fl -P344) = 0.5 (i = 1,...,6) (see Fig. 5. a, 
Fig. 6.c). 
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An=E, 




I(xJ 




Figure 5: The site of the horoball B 3 of the basic horoball arrangements B l 3U 
(i = 1, 2, 3) in the octahedron P344. 



2. Two "larger horoballs" with centers at E 3 and E 5 are tangent at the center 
of the octahedron, while horoballs at the remaining six vertices touch both 
"larger" horoballs. The point of tangency of the above two horoball types on 
segment I(0)E is denoted by I(xi) where X\ = log(2)/2 is the hyperbolic 
distance between 1(0) and I(xi). In this case V\ := Vol(Bi fl P344) = 1 
(i = 3,5) and V 2 := Vol{B t n P 344 ) = 0.25 (z = 0, 1,2,4) (see Fig. 5.b, 
Fig. 6.b). 

3. One horoball of the "maximally large" type centered at E 3 . The large 
horoball is tangent to all non-neighboring sides of the octahedron and it 
determines the other five horoballs touching the "large horoball". The point 
of tangency of the two horoballs along segment I(xi)E is denoted by I(x2) 
where x 2 = — log(2)/2 is the hyperbolic distance between 1(0) and I(x 2 ). 
In this case V 3 := Vol(B 3 n P 344 ) = 0.25, V 5 := Vol(B 5 n P 344 ) = 0.0625 
and Vi := Vol(Bi D P344) = 0.03125 (i = 0, 1, 2, 4) (see Fig. 5.c, Fig. 6.a). 



Due to symmetry considerations it is sufficient to consider the cases when one 
horoball extends from the "midpoint" of an edge until it touches the opposite side 
of the cell. Here we give an analogous argument as in the proof of Proposition 4.4. 
Assume that none of the horoballs covers a "midpoint". Then the packing density 
may be improved until at least one ball reaches a midpoint. Moreover, consider 
that point I(x) is on edge A A 3 = E E 3 . This point is where the horoballs 
Bi(x), (i = 0, 3) are tangent at point I(x) G A I(0). Then x is the hyperbolic 
distance between 1(0) and I(x), and it is analogous to the previous section. It 
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is easy to see that we have to study two different cases to determine the optimal 
horoball arrangement: 

1. x e [0, xi], horoballs P 3 and P 5 touch horoballs Bi (i = 1, 2, 3, 4). 

2. x G [xi, X2], horoball £? 3 touches horoballs B; L [i = 1, 2, 3, 4, 5). 

In the first case the function V(x) can be computed by the following formula 

V(x) := 4 • Vol(B (x) n P 344 ) + 2 • Vo/(-B 3 (a;) n P 344 ) x G [0, x x ]. 
Similarly to the Lemma 4.3, we can prove the following Lemma: 
Lemma 4.9 

V(x) := 4 • Vol(B (x) n P 344 ) + 2 • \/o/(P 3 (x) n P 344 ) = 
= \/ (2-e 2:E + 4-e- 2:E ), ig[0,Ii], 

and the maxima of function V(x) are realized in points 1(0) and I(xi). 

In the second case, similarly to Lemma 4.6, the volume function V(x) is given 
by the following formula: 

Lemma 4.10 

V(x) := 4 • V ie - 2{X - Xl) + V 2 ■ (e 2{x ~ xi) + e - 2(x - xi) ), x € [x u x 2 ], 
and the maximum of function V(x) is achieved at points I(xi) and I(x 2 ). 

By Definition 4.2 and Remark 4.1, as well as by Lemmas 4.5-4.6, it follows 
that the densest horoball packing is realized in three distinct primary horoball ar- 
rangements £> 344 (i = 1, 2, 3). These optimal horoball packings belong to horoball 
parameters s = —1/3, s = 0, s = 1/3 (see Section (3.1)) and are illustrated in 
Fig. 6.a-b-c. 

The maximal density can be computed by the method described in Section 4.2. 
Thus we have proven 

Theorem 4.11 Three different optimally dense horoball arrangements i3 344 , (i = 
1, 2, 3) exist for the octahedral Coxeter tiling (3, 4, 4), which share the optimum 

density 5(B l 3U ) w 0.818808. 
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a. s = -1/3 b. s = c. s = 1/3 

Figure 6: Optimal horoball packings of Coxeter honeycomb (3, 4, 4). 



4.4 The (4, 3, 6) Cubic Tiling 

The optimal packing densities for the cubic Coxeter tiling (4,3,6) can be obtained 
by similar approaches as in the previous two sections for (3, 3, 6) and (3, 4, 4). 

We consider horoball packings with centers located at the ideal vertices of the 
cube honeycomb (4, 3, 6). 

Analogous to the above cases we introduce a projective coordinate system, by 
an orthogonal vector basis with signature (—1,1,1,1), with the following coordi- 
nates of the vertices of the infinite regular cube (see Fig. 7), in the Cayley-Klein 
ball model: 



y/2 V2 1 



\/2 y/2 I, 



V2 1 



£ (l,-^=,^-,-), ^(1,-^- _^_ ), e 2 (1,0,2^-,-- 



v/3' 3 '3 



v/3' 3 ' 3' 

V2 V2 1, 



£ 3 (1,0,0,1), £4(1, T , ■ ;:)/ . 

Using Lobachevsky volume formula for orthoschemes, we obtain the volume of 
one cubic cell HSz05-2H : Vol(P 436 ) = 48 • Vol(0 {4j3t6) ) « 0.507471. Applying 
the definition of the packing density for the case of tiling (4, 3, 6), we obtain: 



S(B i36 ) 



ELiVoi(B t nP i36 ) 



Vol(P, 



(4.7) 



436 J 



where Bi fl P 436 (i = 1, . . . , 8) denote the 8 horoball sectors, one in each vertex 
of the cube P 436 and we assume that the horoballs B L form a packing in H 3 . We 



consider the following four main horoball configurations B 



1.36 



1,2,3,4): 
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>E 7 



Figure 7: The site of the horoball £> 3 of the base horoball arrangements B\ m 
(i — 1, 2, 3, . . . , 8) in the cubic tiling P436. 



1. All 8 horoballs are of the same type and the adjacent horoballs touch each 
other at the "midpoints" of each edge (see Fig. 8. a). The density of this 
packing: 5{B\ m ) « 0.682621. 

2. Two "larger horoballs" with centers at E 3 and E 7 are tangent at the center 
of the cube, while horoballs at the the remaining six vertices touch both 
"larger" horoballs (see Fig. 8.b). The density of this packing: 5(B\ m ) « 
0.682621. 

3. One large horoballs in same types touch each other in the "midpoints" of 
the "diagonals" of cube faces. The other four horoballs touch the adjacent 
"large horoballs" (see Fig. 9. a). The density of this packing: 5(B\ m ) ps 
0.853276. 



4. One horoball of the "maximally large" type centered at E 3 . The large 
horoball is tangent to all non-neighboring sides of the cube and it deter- 
mines the other five horoballs touching the "large horoball" (see Fig. 9.b). 



The density of this packing: 5(B. 



1.--S6 J 



0.853276. 



By using Lemma 4.6 and its corollaries, similarly to the above cases we may again 
write the volume function of the horoball pieces to prove the optimality of two of 
the four limiting cases. Finally, we obtain the following 



Packings for Fully Asymptotic Coxeter Tilings.. 



21 





a. b. 

Figure 8: Locally optimal packings of Coxeter honeycomb (4, 3, 6) 





a. b. 

Figure 9: Optimal horoball packings of Coxeter honeycomb (4, 3, 6). 

Theorem 4.12 Two different optimally dense horoball arrangements B\ m , (i = 
3, 4) exist for the cubic Coxeter tiling (4, 3, 6), which share the optimum density 

5(Bi 36 ) « 0.85327609. 



5 The (5, 3, 6) Dodecahedral Tiling 

The optimal packing densities for the dodecahedral Coxeter tiling (5, 3, 6) is ob- 
tained through a similar approach as in the previous three section, hence we omit 
the details of the derivation. 

We consider horoball packings with centers of horoballs located at ideal ver- 
tices of the dodecahedral honeycomb (5, 3, 6). 
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Analogous to the previous cases we introduce the projective coordinate sys- 
tem, by an orthogonal vector basis with signature (—1, 1, 1, 1), with the following 
coordinates of the vertices of the infinite regular dodecahedron, in the Cayley- 
Klein ball model. The dodecahedron contains a cubic sub-lattice with coordinates 
adopted from the previous section, as well as 12 other vertices obtainable through 
rotation about these vertices. 

Using Lobachevsky volume formula for orthoschemes, we obtain the volume 
of one dodecahedron HSz05-2H : Vol(P 536 ) = 120 ■ Vol(0 (5A(i) ) « 20.580199 .... 
Applying the definition of the packing density for the case of tiling (5, 3, 6), we 
obtain: 

S(Bat) = g^^M , (4.7) 

where B, t fl P536 (i = 1, . . . , 20) denote the 20 horoball sectors, one in each 
vertex of the dodecahedron P536. By using that the dodecahedral tiling has a 
cubic sublattice, we consider the following five main horoball configurations B l 536 

(i = l,...,5): 

1. All 20 horoballs are of the same type and adjacent horoballs are tangent 
at the "midpoints" of each connecting edge. The density of this packing: 
tfOSke) « 0.550841. 

2. Two types of horoballs occur in this packing confguration. Eight larger 
horoballs are centered at the lattice points of the dodecahedra making up 
a cubic sublattice as in packing B\ m . Twelve smaller types of horoballs 
are located at the remaining 12 lattice points. The density of this packing: 
S(Bl m ) » 0.70309. 

3. This packing configuration contains horoballs of 4 types. The cubic sub- 
lattice within the dodecahedral lattice has the same ball configuration as 
£>| 36 , and the two types of the balls on the cubic lattice points uniquely de- 
termine two types of neighboring horoballs. The density of this packing: 
S(Bl m ) « 0.78725. 

4. There are three types of horoballs in this packing. The cubic sublattice 
within the dodecahedral lattice has the same ball configuration as £>f 36 , 
and the larger horoball uniquely determines the type of the remaining 12 
horoballs. The density of this packing: 6(B% 36 ) « 0.78481. 
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5. This limiting case is an extension of packing <Bf 36 . We inflate the horoball 
located at (1, 0, 0, 1) until it touches the non-adjacent side. This packing 
consists of 6 horoball types, 4 of which are on the cubic sublattice and are a 
non-limiting case of the cubic tiling. These 4 horoballs uniquely determine 
two horoball types of the remaining 12 vertices of the dodecahedral lattice. 
The density of this packing: <5(Sf 36 ) « 0.71246. 




Figure 10: Two tangent horoballs in a packing of the (5, 3, 6) honeycomb. 



By using Lemma 4.6 and its corollaries, similarly to the above three cases we 
may again write the volume function of the horoball pieces to prove the optimality 
of two of the four limiting cases, leading to the following theorem: 



Theorem 5.1 The optimally dense horoball arrangement B\ m for the dodecahe- 



dral Coxeter tiling (5, 3, 6) has optimal density 8{B\ 



536 ) 



0.787251. 



6 Conclusion 



Locally optimal horoball packings were studied in this paper, in which different 
types of horoballs were placed at the lattice points of fully asymptotic Coxeter 
honeycombs. We proved that the value obtained by Boroczky and Florian as the 
universal density upper bound for all congruent ball packings in hyperbolic 3- 
space remains the upper bound even if horoballs of different types are considered. 
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Moreover, there are two distinct optimally dense packings for the (3, 3, 6) Coxeter 
tiling. We again encounter the Boroczky - Florian density upper bound as upper 
limits when varying the types of horoballs at the lattice points of the (4, 3, 6) cu- 
bic Coxeter tilings, with two distinct realizations of optimality. For the (3, 4, 4) 
octahedral Coxeter tiling, the optimal packing density is less than the maximal 
value, and there are three distinct configurations of balls yielding the same opti- 
mal value. The case if the (5, 3, 6) dodecahedral tiling we have obtained some 
interesting horoball arrangements with less densities. Table 2 contains a summary 
of the optimal packing densities under our constraints. 
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Table 2. Optimal packing densities for the four fully asymptotic Coxeter tilings 



Schlafli symbol 


Optimal density 


(3,3,6) 
(3,4,4) 
(4,3,6) 
(5,3,6) 


0.853276* 
0.818808 
0.853276* 
0.787251 



*These values are identical to the Boroczky and Florian limit. 

In the future it will be interesting to investigate tilings given various uniform 
conditions on the configuration of the balls in H 3 as well as higher dimensional 
hyperbolic spaces. These studies may show the existence of multiple optimal con- 
figurations for given tilings, similarly as we have observed in H 3 . To the knowl- 
edge of the authors the solution of the above problem is still open. 

Optimal sphere packings in other homogeneous Thurston geometries repre- 
sent another huge class of open mathematical problems. For these non-Euclidean 
geometries only very few results are known HSz07-2L HSzlO-lH . Detailed studies 
are the objective of ongoing research. 
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